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Abstract 

We consider the spherically symmetric metric with a comoving perfect fluid and non-zero 
pressure — the LemaTtre metric — and present it in the form of a calculational algorithm. 
We use it to review the definition of mass, and to look at the apparent horizon relations on the 
observer's past null cone. We show that the introduction of pressure makes it difficult to separate 
the mass from other physical parameters in an invariant way. Under the usual mass definition, 
the apparent horizon relation, that relates the diameter distance to the cosmic mass, remains the 
same as in the LemaTtre- Tol man case. 

1 Introduction 

The modern theory of cosmology is totally based on the general theory of relativity. Underlying nearly 
all modern cosmological is the assumption of a homogeneous FLRW model. Successful as it is, there 
are many unanswered questions about the nature of the matter content of the cosmos, as well as 
the relationship between observations of many discrete sources, which are necessarily averaged in 
some sense, and general relativity, which assumes the metric and the matter are described by smooth 
functions. A significant complication arises from the fact that observers are looking down earth's past 
null cone, so we don't know the state of the universe at any one time or the history of any one worldline. 

When attempting to verify the question of homogeneity — something we should do in the near 
future — there are two distinct aspects. Isotropy about us is relatively easy to check — either the 
observations (at a given redshift) are (statistically) the same in different directions or they are not. 
One doesn't need to know the spacetime evolution to check this. But establishing radial homogeneity 
requires that we interpret observations in terms of a model that is not already assumed to be homo- 
geneous. The observations we make are affected by the bulk equation of state (which fixes the cosmic 
evolution), and also by any radial inhomogeneity that might be present. These affect the shape of the 
past null cone and hence the time at which we observe each worldline. The evolution of the sources 
properties must also be known, and theories of source evolution that have been validated or derived 
assuming a homogeneous model cannot be used to verify homogeneity. 

*ALFALN001@uct.ac.za 
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The "observational cosmology" series and related papers [10, 29, 30, 31, 32, 21, 22, 4, 2, 3, 28, 
1, 14, 5] considered how the spacetime geometry could be determined from cosmological observations, 
starting with the classic by Kristian and Sachs [18], and using key concepts in [11, 33]. An important 
construction is the use of observer coordinates, based on the past null cone. Recent work on obtaining 
the metric of the cosmos from observations [25, 12, 13, 20, 23] has focussed on the practicalities of 
a numerical procedure for converting observations into metric information. To date this investigation 
has assumed spherical symmetry, as the key first step away from homogeneity towards the far more 
complex general case. It was pointed out [13] that the apparent horizon — the locus where the diameter 
distance is maximum — has important observational consequences, that enable us to determine the 
cosmic mass at that radius. This is due to a simple relationship between the maximum in the diameter 
distance, the cosmological constant, and the mass within a sphere of that size, that is unique to that 
distance. This relationship allows us to check for systematic errors in the observational data, and 
partially correct them [23]. 

As early as 1933, LemaTtre [19] considered the general, diagonal, spherically symmetric metric, 
imposing a comoving, diagonal, matter tensor, but allowing the radial and tangential pressures to be 
different, and the cosmological constant to be non-zero. He obtained the conservation equations, and 
from them obtained an expression for the mass. He reduced the field equations to a system of first 
order differential equations (DEs) constraining the initial configuration and fixing the evolution. He 
then went on to consider a large variety of very interesting cases and applications — see [17]. This 
paper was well ahead of its time, and is little appreciated even today.^ 

In 1964, Podurets [27] wrote out the field equations for a spherically symmetric "star" consisting of 
a comoving perfect fluid (with isotropic pressure) and zero cosmological constant. He solved them down 
to a slightly simpler set of evolution and constraint DEs, gave a brief thermodynamic interpretation of 
the equation relating pressure to the time derivative of mass, and presented only a sketch of how a 
numerical procedure would work. In order to gain a basic interpretation of the equations, he showed 
how the Euler, continuity and Poisson equations are obtained in the Newtonian limit. 

That same year, Misner and Sharp [24] considered the same kind of star, pointing out it cannot 

include heat conduction or a flux of radiation. They discussed the thermodynamics in some detail, 
and from this, they obtained the the same evolution and constraint DEs — the relativistic Euler and 
continuity equations. They defined a particle number and showed it is conserved. They imposed an 
origin condition at zero radius, and a comoving boundary condition where the metric joins to a vacuum 
exterior. 

Cahill and McVittie's [8] derivation of the spacetime is essentially the same as the above two, 
except they don't assume a diagonal matter tensor. They provide several justifications for the choice 
of mass function, which we discuss below. They go on to consider the perfect fluid, 'negative mass 
shells',^ and some specific metrics. 

See [16] for an excellent review of this and many other inhomogeneous cosmological models. 

Given the significance of the apparent horizon result, it is important to check how general the 

^Possibly the size of the paper was a hindrance to it's appreciation — it contains enough material for several 

important papers. Another aspect that possibly makes it difficult to penetrate at the start is the rapid introduction 
of a lot of new variables, the as and /3s. The initial diagonal metric seems to have signature [+ + ++], which means 
that the spatial metric components, af, a| & 03 are really negative (in the signature adopted later) so ai to 03 are 
imaginary. This makes the definition for the mass m = — 47rz4> rather peculiar, because $ is imaginary, and it also 

introduces unexpected signs in <l>'s definition (2.12). Nevertheless, from section 3 onwards, the signature [ h] 

and the notation are quite normal and the solution equations, e.g. (3.5), take an easily recognisable form. 

^Since they assume p Sz p are positive, these shells would appear to be the region beyond a maximum of R in 
the spatial sections of constant t, where both R and M are decreasing towards a second origin. 
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relationship is. The main content of this paper can be express as follows: In section two we will 
present the general LemaTtre metric, review how the solution of Einstein's field equation is reduced 
to a system of differential equations (DEs), and then show how one can construct a LemaTtre model 
using appropriate numerical integration of the DEs, and finally we will show how the LTB and RW 
special cases are obtained. Section three will focus on how the definition of mass is affected by the 
introduction of non zero pressure and cosmological constant. In section four we will extend the apparent 
horizon relationship to the case of non-zero pressure, and discuss the past null cone with its observable 
quantities. We end this paper with the conclusions. 



2 LemaTtre Metric 

We consider an inhomogeneous, spherically symmetric spacetime, filled with a perfect fluid, in which 
the coordinates — {t, r, 9, 0) are comoving with the matter flow. The metric may be written as 

ds^ ^ -e^"" dt^ + dr^ + R'' dn\ (1) 

where a = a{t,r), A = A(t, r) are functions to be determined, R = R{t,r) corresponds to the areal 
radius, and dfl^ = dO^ + siv? d(j? is the metric of the unit 2-sphere. The energy momentum tensor 
is given by: 

Ti"' ^{p + p) vy + g^^^v , (2) 

where p — p{t,r) is the mass-energy density of the perfect fluid, p — p{t,r) is the matter pressure, 
and u'^ — {e~^, 0, 0, 0) is the fluid four-velocity. This is less general than the full LemaTtre metric as 
we assume isotropic pressure. 



2.1 Field Equations 

The EFEs, G^"^ = kT^^" — gf^^A can be reduced to the following set of equations: 



e^G'^^(^-^a'-^X \ e-^'^^O. (4) 
^ R R R ' ^ ^ 

,A r^rr [ ^'^ , A _x ( 2R R^ 2i? . \ 1 



+ (^a-^--\ + -\a-^\-->?\ e-'"' ^Kp-K. (6) 
yR R 2 2 2R A ^ ^ ^ ' 

Where the dot means a derivative with respect to t, and prime means a derivative with respect to r. 
We use geometric units, G = 1 = c, so that k — Stt. The conservation equations V^T'^'^ — are: 

V,T*^^A+^ + ^=0 (7) 



(p + p) (p + p) R 
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V^T^'^ = a' + = . 



(8) 



To solve these equations, we multiply Eq (3) by E?R' and use Eq (4) to eliminate the term that 
contains A, which produces 



d_ 
dr 



R + RR'^e-^" - RRl'^e'^ - -AR^ 

3 



KpR^R! 



(9) 



Multiplying Eq (5) by R^R and using Eq (4) to eliminate the term that contains a', shows Eq (5) can 
be rewritten as 



d_ 

di 



R + RR^e'^" - RR'^e-^ - -AR^ 



-KpR^R 



(10) 



The term in square brackets is related to the total mass-energy of the system, M, interior to a comoving 
shell of constant r, and is normally defined by 



^ ^ -2a _ -A ^ 1 _ 1^^2 

R 3 



(11) 



The justifications for this will be discussed in section (3). With this definition, Eqs (9) and (10) can 
be rewritten as: 



2M' 

Wr' ' 

2M 
~R^ 



(12) 
(13) 



These two relationships will be discussed in §3. Eq (11) may be rearranged as an evolution equation 
for the model 



i? = ±e^/ — + / + 



R 



3 ' 



(14) 



where 



f{t,r)^R%-'-l, 



(15) 



acts as the curvature term, or twice the total energy of the particles at r (analogous to / in the 
LT model). However, the solution for R{t,r) cannot be directly obtained from this, because of the 
unknown functions A, a and M. 

The metric variables gu and grr can be obtained by integrating Eqs (8) and (7) as follows: 



const t const t 

r^..„_ /',(*^ dp. (16) 

J ro 



(P + P) 



Po 



{P + P{p)) 



and 



const r 
P 



A = An(r 



''Po 



{p + p{p)) \Ro 



(17) 
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where ao{t) and Ao(r) are arbitrary functions of integration. Typically, we would choose tq = to be 
the origin, where R{t,0) — 0. Having solved (3)-(5) and (7)-(8), the 9-9 field equation (6) is now 
satisfied, since 

AR' SRl 2) \Ar) 



AR' 8i?y 2 / V ^R' 



and to get the right hand side of (6) we must also add 



= ^ V.T- - ^ V.T- . (19) 



2.2 Constructing a LemaTtre Model 

To actually generate a LemaTtre model in the general case requires a numerical integration of the DEs. 
In order to do that, we first need to specify the free functions, and then integrate the DEs along 
constant t or r paths. Below is a summary of the procedure. 

• Choose an initial time slice Iq, and make a gauge choice R{to,r) = Ro{r), for example Rq = r. 
This means Rq is also known. 

• Specify po{r) = p{to,r) on that initial time slice. 

• Select a specific equation of state, p — p{p), thus giving po — p{to: f), on the initial surface. 

• Integrate Eq (12) along constant t — t^ 

const t 

Mo= r ^!^5^dr, (20) 



to calculate Mq = M{tQ^r) everywhere on the initial time slice. Alternatively, one may specify 
Mo(r) and determine PQ{r) from it. 

• Choose A(to,''), thereby fixing the free function Ao(r). Note that by (15) this is equivalent to 
choosing the geometry /o(r) on the initial surface, so Ao(-R) = 0, eg = 1 is not necessarily a good 
choice. One may instead prefer to choose /o(r) and then calculate Ao(r) from = it!o /(l+/o). 

• Choose cr(t, ro) along the central worldline Tq = 0, which we can take to be identically cro(t). 
This relates the time coordinate to the central observer's proper time. 

• Integrate (16) along t = to to get a{t(i,r). 
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We have now determined all the functions on the initial time surface, which means we can obtain their 
r derivatives too. We next wish to evolve the metric forwards in time, along the worldlines of constant 
r. This can be done as follows: 



• Equation (14) gives R everywhere on our initial time slice. 

• The M DE can be obtained from Eq (13) as 

M = i-- . 



(21) 



Eliminating a' between Eqs (8) and (4), and then substituting for A from (7), produces the 
following DE for p 



R 



R' 2R 



(22) 



Having chosen the equation of state^ p — p{p), an equation for p follows, 

dp . 
dp 



(23) 



Eq (7) provides a DE for A, which may be used as is, or combined with (22) to give 



p'R 
R' \{p + p) 



+ R' 



(24) 



• Having these 5 coupled DEs, and the initial values ont = to, we can create a numerical procedure 
to integrate the DEs in parallel, thus giving R{t,r), M{t,r), p{t,r), p{t,r) and \{t,r) every- 
where. Eq (12) can be used as a cross-check on the results. Note that the spatial derivatives R! , 
p' and M' will be needed at each step. 

• Finally, cr{t,r) is obtained from Eq (16) by integrating along each slice of constant t. 



Notice that we have chosen 4 functions, Ro{r), po{r), Xo{r) and ao{t), as well as the equation of state 

P{p)- 



2.3 Special Cases 

The LemaTtre metric contains the LemaTtre-Tolman and FLRW metrics as special cases. 

The spherically symmetric inhomogeneous dust cosmology of LemaTtre-Tolman is obtained when 
the pressure is zero. Setting p = \n (8) gives a = ao{t), and we are free to set a = 1 so that t 
becomes the comoving proper time, i.e. the natural choice of time coordinate. Eq (13) gives 

M - ^ M = M(r) = Mo(r) . (25) 
^More general equations of state would require a slightly different procedure 
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Similarly, (7) and (17) simplify to 



A = Ao(r) -21n 



4 In 



e ' = 



P_ 

Po 



R_ 

Rq 



R' 



;i + /o) 



(26) 
(27) 



where we used (12) and then (15), so clearly / = foir) — f{r) too. Therefore the evolution equation 
reduces to the familiar 

• 2 2M ^ A7?2 

^^^ + f + ^- (28) 

On the other hand, the standard cosmological model, the homogeneous RW metric, is obtained by 
setting to zero the spatial variation of any physical invariants. Again p' — allows us to choose a = 1, 
and requiring the constant t 3-spaces to have a constant curvature form requires gij — S'^{t)gij{r), 
i,j = 1,2,3, with the canonical r coordinate choice giving 



ds^ = -df + S{t) 



1 — kr'^ 



(29) 



3 Definition of Mass 

In this section we will consider the justification for naming M the gravitational mass felt at comoving 
radius r. First we review the arguments given by some earlier authors, especially Cahill and McVittie, 
then we consider a geodesic deviation approach. 

Lema?tre [19] does not really justify calling M the mass, other than an implied comparison with 
Newtonian equations. Podurets [27] (with A = 0) merely says "it is not difficult to see that" (11) 
gives "the total mass of matter contained in the interval (0, r), the mass being defined in terms of the 
gravitational field created on the boundary", and then calls it a definition of mass. Misner and Sharp 
[24] (also using A = 0) write (12) in the form 

M = / p(l + ife-^'^ - ^) d V , (30) 

Jv(r) V R J 

pointing out that this M contains contributions from the kinetic energy and the gravitational potential 
energy, as well as the matter density p. (We note that the contributions f/^ = R^e~'^"' and 2M/R 
are twice the two energy values. We also note that the term M appears in both sides, which causes a 
difficulty in interpreting this equation.) 

Cahill &i McVittie [8], assuming zero A, define the mass in terms of a Riemann component 

McMit, r) = V = I + |! - ^ j . (31) 

which involves only first metric derivatives, and is independent of any rescaling of the t r coordinates. 
Their justifications are: (i) if joined to a Schwarzschild exterior, M must equal the exterior mass; (ii) 
the Bianchi identities lead to generalisations of (12) &i (13), viz 

2M' = kR^ {Tl R! - R) , 2M= kR^ (T; R - T[ R!) , (32) 
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but the roles of the various terms in the interpretation are not discussed; (iii) the proper acceleration 
of i? on a comoving worldline is 

u^V^iu^V.R) = -RT; - — + R'e-^X' (33) 

which contains a Newtonian-like gravitational force term —M/R^, as well as the pressure acting on 
the comoving shell, and a term that's hard to interpret; (iv) they show that the mass flow vector 

r = -^^^(M', -M, 0, 0) , r,^ = , (34) 

is conserved.'* (For the case of the Reissner-Nordstrom metric, with mass and charge parameters m 
and q, they find that 

M^m- = m - 27r / R^ dR , (35) 

R Jr 

where E is the electric field. The two masses agree at i? = oo, but the effective gravitational mass M 
decreases as the charge is approached, because the energy density in the E field that is outside radius 
it! is not included in M. As is well known, this energy density diverges close to the charge, so M can 
go negative.) 

Aspects of the above arguments are clear and easy to follow, and for the A = = p case they 
are convincing. It is not quite so easy to be certain we've properly understood the various relativistic 
corrections when p ^ 0. A strong argument is that (12) is the same as in the LT model, and it applies 
at all times. In other words, on every constant time 3-space, the gravitational mass is the integral of 
the density with a curvature (or gravity) correction, d^V being the proper 3-volume element: 



M 



[ p^/l + fd'V. (36) 

Jvir) 



A similar treatment of (13) might seem to involve more complicated contributions from the gravity 
(curvature). If we consider all the matter within the comoving sphere of constant r, then the work 
done by the pressure on the matter outside the boundary as it expands is dW — pAdL, where A is 
the area of the boundary and dL its displacement. Now if, following traditional Newtonian thinking, 
dL is taken to be the physical distance the boundary has moved, then this is 

dW =pA^( [ (37) 

dr ) dt (38) 



dt \Jo 

^) (4rf)Afr_J= 

KRmj ^ ^dt v7o \/r+7 

R' 



_ _M d f r 



VT+f 



dr) dt . (39) 



which differs from the loss of mass-energy in the interior. If instead we consider how much the boundary 
has expanded, not obviously the correct thing to do, since it has more to do with increase in surface 
area than radial distance moved, it gives the familiar Newtonian result:^ dL — Rdt leads to 

dW = -Mdt. (40) 



'*Their (3.5) and (4.3) also seem to use a comma for a covariant derivative. 

^This thinking does not extend to (36), which is unavoidably an integral over a volume, whereas (40) is really 
about the boundary only. 



8 



Another problem arises once the cosmological constant is introduced, because Cahill &i McVittie's 
definition (31) then involves more than just the mass, 

|i?V = M+l^. (41) 

Consequently, in the following we investigate whether new invariant expressions can provide definitions 
of M, A and other quantities independently. 



3.1 Geodesic Deviation Equation 

We here look at another method that might clarify whether M is the total mass-energy, based on the 
geodesic deviation equation (GDE). The GDE measures the relative acceleration between two nearby 
geodesies, which is itself a measure of the tidal effects of gravity. For a congruence of geodesies, 

^-^^-R^^^^U^eU^ (42) 

where R^upx is the Riemann tensor, Lf^ is the geodesic tangent vector dx^/dr, r is proper distance or 
time along the geodesies, and is the geodesic deviation vector. Contracting this equation with 
produces the scalar 

A(e, -iJ^- -R^r^px e e , (43) 

If has unit magnitude, then A is the size of the relative acceleration of the chosen geodesies, and 
hence it is a measure of the strength of the tidal effects. Now the Riemann tensor determines the tidal 
effects due to local as well as distant matter. If one wishes to remove the effect of local matter, then 
one may consider the quantities 

-C^upxU''eu\ (44) 

and 

Bie, = -C,,,, e e , (45) 

where C^^px is the Weyl tensor. Since we only wish to "feel" the tides in a small neighbourhood, and 
we don't need to propagate the geodesies, we are free to choose U'^ and at will.^ 

In order to make these measures meaningful, we construct a set of canonical vectors using defini- 
tions that are invariant for any spherical metric. In spherical symmetry, the "areal radius" or "curvature 
coordinate" is the metric component that multiplies the unit 2-sphere, dJl^. In (1) it is R. Then the 
unit timelike vector that follows constant R is, 

uf^u^ = -1 , u^V^R = ^ 

-fr-A/2 _ 

— (r',-R,0,0) ■ (46) 



Using this, we define a canonical unit spacelike vector v'^ in the radial direction that is orthogonal to 
u'^ and satisfies 

v^'v^, - 1 , v^M/, = = = v"^ 



^More general measures of the curvature may be defined via the parallel transport equation, i.e. using the 
quantities Wa R°'bcd Z'^, but for the spacetimes under consideration, all the Riemann components with more 

than 2 different indices are zero. 
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-i?e-2'", R'e-^, 0, 



(47) 



It is evident that the spacelike or timelike character of u'^ &i v" will flip if (e~'^i?'^ — e^'^ R^) changes 
sign. We then define unit vectors in the 9 and directions that are also orthogonal to and , and 
tangent to the constant R spheres, 



w'^w,, = 1 



u^Wn = v'^w^ = , 



(48) 



and 



z^'V^R = 



u^z,^ 



0, 



0,0,0, 



Lastly, the radial incoming null vector k'^ is given by 

k^'k,, ^ ^ ^ k't' k^" ^k {6^/^,-6", 0,0) , 



(49) 



(50) 



where k is undetermined by the above conditions. In addition, one could use the unit eigenvectors of 
the Einstein tensor (G^^ - %^) = which are e~^5f , e~^/2(5^, R~^5^, (Rsme)-^5'^, but it turns 
out these do not extend the range of independent scalars we can define. 

The scalars A &i B, when combined with any canonical vectors, are determined solely by the 
metric, and therefore represent a property of the spacetime. Below we use the notation 



AH = e-^R'^ - e^^R? = 1 



2M KR^ 



R 3 ' 

AHm = e-^''^R! - e^R , AHp = e~^/^R' + e^R . 



Evaluating the scalars A and B using the various combinations of the above vectors gives: 







A{u'',v'') 






B{u^',v'') 


Aiu'', 




= Aiu'', z") 


Biuf" 


wn 


= B{u>',z') 


Aiu" 


X) 


= Aivf", r ) 


B{u^ 


,r) 


= fi(^;'^,r) 


Aivf" 




= Aiv'', z") 


Biv^" 




= B{v^,z'') 



2M A K. 
i?3+3 


+ P) 




2M Kp 






i?3 3 






up k{p + p)R' 




M 


2 2{AH) 


~ R^ 


Kp M 






Y~R^ 








'2M 


A 


(AHp) 




k^e^+^^iAHm) 


"2M 


up 


(AHp) 


i?3 





K 
2' 



Kp K{p + p)R'^e~^ M A 
'2 2{AH) ^R^ ~ 3" 

M Kp 



(51) 
(52) 

(53) 
(54) 
(55) 
(56) 
(57) 

(58) 
(59) 
(60) 
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^K.^T--^-t (61) 

= (62) 

E{w>' .k") = B^z'' X) = ^ (64) 

We notice that all the Bs are multiples of — 2M/ B?, while the are more diverse. Combining 
the above equations, we can define the following 

^K,»^)-^K,™T = ^-|^-^ (65) 

2M A 

-^K, ^'^) = ^ + 3 (66) 

A{v^, w") - ^(m'', w") - ^(m^, = - A (67) 

-2\a{uI', v"") + ^(w'^, z'') } - /«(p + p) (68) 

^{^(^;^^^;'^) -^««;^) -^(w^i;^)} = ^ + y (69) 

-^(«^, ^;'^) - 2A{w^, z") - Aivf", w") + ^(m'^, w") = np + K (70) 
^ {^(^^ ^^) - ^(^^^ w^) + ^(i.^ t;'^) - A{w^, z^) } = z'^) = ^ - ^ . (71) 

However, there are no combinations that give the value of just M/R^, or np, or np or A. In the case of 
zero A, though, Eqs (66), (67), (70) do give M/R^, np and np. Similarly, in the case of zero pressure, 
Eqs (67), (68), (69) do give A, Kp and M/R^. This suggests that there is no definition of mass based 
on physical invariants in the general case. 

It appears from the results above, that we cannot separate the mass, the cosmological constant, 
the density and the pressure from each other, and so we cannot create a unique definition of mass 
based on geometric invariants of the metric in the general case. In contrast, if p = or A = (the 
LT and Misner-Sharp-Podurets cases), then the remaining quantities are uniquely defined from these 
scalars. This has important implications. The results of Cahill and McVittie do not generalise to the 
case of non-zero A, and the justification for calling the M of (11) the gravitational mass is weaker 
than was thought. It is based solely on comparing (12) and (13) with a few special cases. 



4 Past Null Cone, Apparent Horizon, and Cosmic Mass 

We consider the past null cone (PNC) of the observation event t = to, r = ro = 0. For any quantity 
Q{t,r), its value on our PNC will be indicated with a hat: Q — Q{r) — Q{i{r),r), or for expressions 
a square bracket with subscript "A" will be used. For the metric (1), the path of an incoming radial 
light ray is given by 

f- = , (72) 

and we define the solution that reaches {to,ro) to be t = i{r). 
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The redshift of comoving sources on that null cone is given by the ratio of the light oscillation 
periods T measured at the observer, o, and the emitter, e. 



(73) 



For the two successive wavefronts passing through the events B and A or D and C respectively on a 
worldline of constant r, and a neighbouring one at r + dr, the change in the light oscillation period T 
over a distance dr is given by 



dT 



dr 



dr — 



c 



dt 
dr 



dr , 



(74) 



The first term on the right of Eq (74) can be written as a Taylor expansion about A, so that 



dT 
dT 



dt 
dr 
d_ /dt 
dt Vdr 



+ 



d fdt 



dt Vdr 



dr 



T + 



d /dt 



dr V dr 



(0) } dr - 



dt 
dr 



dr 



Therefore Eq (74) can be integrated down the PNC to give 



ln( ^ 

T. 



^(^)'^' 



so that, for the LemaTtre metric, Eqs (72) and (73), give the redshift as 



ln(l + z] 



rr. Q /gA/2 



dr 



(75) 
(76) 



(77) 



(78) 



With any cosmological sources there are primary observable quantities such as redshift z, angular 
diameter 6, apparent luminosity i, and number density in redshift space n. Associated with each of 5, 
£, and n is a source property, true diameter D, absolute luminosity L, and mass per source /x, which 
may evolve with time, and therefore vary with redshift. These combine to give the luminosity distance 
(11, the diameter distance dn and the mass density in redshift space, iin. 

We define the diameter distance of a source as the ratio between the true diameter D and the 
angular diameter 6. Consequently, in a spherical metric, it corresponds to the areal radius evaluated 
on the PNC, i.e. 



D 
J 



dn^R 



(79) 



Given the absolute luminosity of a source L, and the apparent luminosity I (or m and m, the apparent 
and absolute magnitude), then the luminosity distance is 



(80) 



where diQ is 10 parsecs. According the the reciprocity theorem [11, 26, 9], the luminosity distance di 
is related to the diameter distance do by 



do = dL{l + zf . 



(81) 
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Working in redshift space, {z,9,(f)), let n{z) be the density of sources, that is the number per 
steradian per unit redshift interval. Suppose that there are dN sources in solid angle du — smOdOdcf) 
between redshift z and z + dz, and that ijl{z) is the mean mass per source, then 



dM. — n dN — fj,ndujdz 



(82) 



The 3-d volume d^V of proper space that encloses these sources at the time of emission t^, as measured 
by comoving observers u'^, is given by d^V = e^^"^ dio dr , so that the mass in the fluid element is 



dM^ pe^/^R^dudr . 



Therefore the relationship between n and p is given by 



pn 



dr 



(83) 



(84) 



The apparent horizon is the locus where the observer's PNC reaches it's maximum areal radius 
(diameter distance), which we denote z — z^, R — Rm, and obviously we have [d/2/dr]^ = 0. The 
total derivative of R along the ray is 



cLR 

dr 



R' + R 



dr 



and using Eqs (14) and (72) to replace R Sj. di/dr shows 



dr 



+ R'^e-^ - IH 



R 



= 



6M„ - mm + ^Ri 







(85) 



(86) 
(87) 



at Rm- This is exactly the same as the result in [13], showing that the non-zero pressure does not 
affect the apparent horizon condition. Note too that M and A appear in exactly the same inseparable 
combination as in (66), so (87) can be written R^ — —1/A{w'^,z^)\m- 



5 Conclusions 

We have considered the Lemartre metric, which describes a spherically symmetric perfect fluid distribu- 
tion in comoving coordinates, and its familiar reduction to a system of DEs. We presented the solution 
in the form of an explicit algorithm for numerical calculation, clearly listing those functions that need 
to be specified on an initial surface, and on a chosen worldline. 

The concept of the total gravitational mass within a given comoving sphere has been reviewed 
in the case of both pressure and cosmological constant being non-zero, and, in addition to the usual 
relativistic corrections whose role is hard to pin down, it was found that A appears in the Cahill 
McVittie definition of M. In seeking for alternative definitions or justifications, we considered a number 
of invariants of a spherical metric, based on the geodesic deviation equation and a set of canonical 
vector fields. We found that the physical variables M, p, p and A always appear in combinations of 
two or more, and it is not possible to separate them using methods based on invariant properties of the 
metric. Since M is not just the integral of the density, but contains a contribution from the curvature. 
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and since the relation between M and the work done by the pressure is not straightforward, we are 
therefore left with some ambiguity in the definition of the effective gravitational mass. 

We looked at the past null cone of a central observer, and how the areal radius (diameter distance) 
varies along it. The apparent horizon relation found in [13], that is the relationship between the 
maximum in the diameter distance Rm, A, and the mass M within that sphere — the "cosmic mass" 
— has been shown to hold good for the LemaTtre model. This relationship, previously only shown 
for the zero pressure case, holds only at the apparent horizon, and is independent of any intervening 
inhomogeneity. It provides a way of measuring the combination (M^ + AR^/6) on gigaparsec scales. 

Now it is known [25] that the Lemartre-Tolman model can fit any given observational functions for 
diameter distance and number counts versus redshift. A generalisation of this theorem to the LemaTtre 
model would seem to offer more flexibility, and thus the possibility of fitting other data. 
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A Barotropic Equation of State 

For the case of the commonly used equation of state, p = wp, with w the equation of the state 
parameter, the LemaTtre model has a more explicit solution. Eq (16) and (17) can be integrated as 
follows 

const t 

^ p{l + w) (l + w) \ptoJ ^ ^ 

where pto = p{t,ro), and 

const r 

A Ao(r) /•* p'dr f R \ -I . f pR^ 



, 21n =- -In ^ , (89) 



where Ror = R{to,r) etc, so that 



„ \ -w/{l+w) / r,T^^ \ -1/(1+"') 
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